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Abstract. The stellar and dust content of spiral galaxies as function of radius has been investigated using near- 
infrared and optical broadband surface photometry of 86 face-on spiral galaxies. Colors of galaxies correlate with 
the azimuthally averaged local surface brightness both within and among galaxies, with the lower surface brightness 
regions being bluer. The colors formed from different passband combinations correlate strongly indicating that they 
probably arise from the same physical process. 

A 3D radiative transfer model was developed to calculate the effect of dust absorption and scattering on the luminosity 
and color profiles of galaxies. Stellar synthesis models were used to investigate the effects of the star formation history 
and the metallicity on the broadband color profiles. Combining all optical and near-infrared data shows that the color 
gradients in this sample of face-on galaxies are best explained by a combined stellar age and metallicity gradient across 
the disk, with the outer regions being on average younger and of lower metallicity. Dust reddening probably plays 
only a minor role, as the dust models cannot produce reddening profiles that are compatible with the observations. 
The observed color differences implicate substantial M / L\ differences, both within galaxies and among galaxies. The 
variations are such that the "missing light" problem derived from rotation fitting becomes even worse. Late-type 
galaxies (T>6) have lower metallicities and are often of younger average age than earlier types and have therefore 
an entirely different M/L\ in most passbands. The near-infrared passbands are recommended for studies where the 
M/Lx ratios should not vary too much. 
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1. Introduction 

For many years broadband colors have been used to ob- 
tain a basic insight into the contents of galaxies. Broad- 
band photometry is relatively easy to obtain and gives an 
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immediate impression of the spectral energy distribution 
(SED) of an object. Broadband colors are particularly ef- 
ficient when used for statistical investigations such as this 
one. Colors have been used to estimate the stellar pop- 



ulatio ns of g alaxies (e.g. S earle et al. [1973y Ti nsley [1980 ; 
Frogel |l985fc Peletier |198§ Silva & Elston |l994|) and it has 



been suggested that colors can give in formation abo ut the 
dust content of galaxies (Evans 1994 ; Peletier et al. 1994 , 
1995 ). In this paper I use radial color profiles to investigate 
the stellar and dust content of galaxies. 

The problem of determining the stellar content of gal- 
axies from integrated SEDs has been approached from 
two sides, often called the empirical and the evolutionary 
approach (for a review, see O'Connell 1987). In the first 
method, stellar SEDs are fitted to the observed galaxy 
SEDs (Pickles [l985fc Peletier ^989)). This method works 
only if one has spectral (line) information. Generally, the 
broadband colors of a galaxy can be explained by a combi- 
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nat ion of the SEDs of tw o or three types of stars (Aaron- 
son 1978 Bershady 1993). In the second, more theoretical 



approach, stellar SEDs are combined, using some knowl- 
edge of initial conditions and evolutionary time scales of 
different stellar populations, to produce evolutionar y stel- 
lar population synthesis models (for reviews Tinsley |l980t 
Renzini & Buzzoni 1986; more recent models are e.g. Buz- 
zoni 



198$ Bruzual & Chariot [1996[ Worthey |199J) 



The papers of Disney et al. ( |1989| ) and Valentijn ( [1990D 
have renewed the debate on whether spiral galaxies are op- 
tically thick or thin. Broadband colors of galaxies can be 
used to examine this problem, because the dependence of 
dust extinction on wavelength causes reddening. This can 
be used to measure extinction at a certain point through 
the disk using a galaxy or anot her object behind it (An- 
dredakis & van der Kruit 1992 ) or to measure extinction 
within a galax y, for instance acros s a sp iral arm dust lane 
(Rix & Rieke |1993| ; Block et al. |l994[ ). To measure the 



global dust properties of a galaxy by reddening one can 
use the color profile. If one assumes that dust is more 
concentrated towards the center (just like the stars), the 
higher extinction in the center produces a color gradient 
that makes galaxies redder inwards (Evans 1994; Byun et 
al. |1994j). 



Most previous studies investigating stellar popula- 
tion and dust properties of galaxies used their integrated 
broadband colors. The use of surface photometry colors 
is less common, as it is easier to compare integrated pho- 
tometry than surface photometry for large samples of gal- 
axies. Integrated photometry samples the bulk properties 
of galaxies, but because the light distribution of galaxies is 
strongly concentrated, one effectively measures the colors 
of the inner regions of galaxies. The half total light radius 
of an exponential disk is ~1.7 scalelengths, while luminos- 
ity profiles are easily traced out to 4-6 scalelengths. There- 
fore, half of the light in integrated colors comes from an 
area that is less than 1/5 of the area commonly observed 
in galaxies (say within £*25)- 

Our knowledge of the star formation history (SFH) 
and the dust content of galaxies improves when we start 
looking at local colors instead of integrated colors. A 
first improvement is obtained by using the radial color 
distribution (i.e. the color profile) of a galaxy. This has 
been commo n practice for ellipti cal g alaxies (e.g. Peletier 
et al. 1990a ; Goudfrooij et al. 1994 ), but not for spi- 
ral galaxies, because elliptical galaxies are assumed to 
have a sim ple SFH and low dust content (but see Goud- 
frooij 1994 ) opposed to spirals. Color studies of spiral gal- 
axies have been concentrated on edge-on systems, in the 
hope to be better able to separ ate th e dust and stellar pop- 
ulation effects (e.g. Just et al. 1996 ). Even more detailed 
information about galaxies can be o btained by t he use 
of azimuthal profiles (Schweizer 1976 ; Wevers et al. 1986 ) 
and color maps, but these techniques require high reso- 
lution, high signal-to-noise observations and are hard to 



become available (Cornell 
al. 1992, Giovanelli et al. 



parameterise to global scales, which means that they can- 
not be used in statistical studies. 

Due to the large variety of galaxies, statistical stud- 
ies of galaxies require large samples. The introduction 
of CCDs into astronomy made it possible to obtain for 
large samples of galaxies accurate optical surface pho- 
tometry in reasonable observing times (Kent 1984). Very 
large data sets of CCD surface photometry have recently 

1987| Han |1992| ; Mathewson et 
1994). Unfortunately, most of 
these samples are observed in only one or two pass- 
bands. Furthermore, the surface photometry is often re- 
duced to integrated magnitudes and isophotal diameters 
to study extinction effects with an inclinatio n test or to 
study the Tully-Fishcr relation (Tully & Fisher 1977 , here- 
after TF- relation). Fast plate measuring machines have 
also produced surface phot ometr y of large sets of galaxies 
(e.g. Lauberts & Valentijn |l989] hereafter ESO-LV), but 
again only in one or two passbands. 

Since near-infrared (near-IR) arrays have become 
available only in the late eighties, near-IR surface photom- 
etry is available for only a few so mewh at larger samples 
of spiral galaxies. Terndrup et al. ( |l994|) observed 43 gal- 
axies in J and K, which w as complemented w ith r pass- 
band photometry of Kent ( |l984j |l986| , |l987|) . They ex- 
plained the observed colors mainly by population synthe- 
sis an d inv oked dust only for the reddest galaxies. Peletier 
et al. 1994 imaged 37 galaxies in the if-passband and com- 
bined their data with the photometry of the ESO-LV cat- 
alog. They explained their surface photometry predomi- 
nantly in terms of dust distributions and concluded that 
spiral galaxies are optically thick in the center in the B 
passband, under the assumption that there are no popu- 
lation gradients across the disk. 

There are two sets of observations that allow a direct 
physical interpretation of color gradients in spiral galaxies: 
1) The current star formation rate (SFR) as measured by 
the Ha flux has a larger scalelength than the underlying 
older stellar population (Ryder & Dopita 1994 ) . There are 
relatively more young stars in the outer regions of spiral 
galaxies than in the central regions. This will be reflected 
in broadband colors of spiral galaxies. 2) From metallic- 
ity measurements of H n regions it is known that there are 
clear metallicity differences in the gas among different gal- 
axies and that there are metallicity gradients as function 
of radius within galaxies (Villa-Costas & Edmunds 1992; 



Zaritsky et al. 1994). If the metallicity gradients in the gas 
are also (partly) present in the stellar components, the ef- 
fects might be observable in the broadband colors. Stellar 
population synthesis models incorporating both age and 
metallicity effects are needed in the comparison with ob- 
servations of radial color gradients. 

Broadband photometry is often assumed to trace bary- 
onic mass, and the transformation from light to mass is 
performed by postulating a mass-to-light ratio (M/L\). 
Both dust extinction and differences in stellar popula- 
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tions will influence M/L\ ratios, most notably in the bluer 
optical passbands. Color differences, among galaxies and 
locally within galaxies, will translate in different M/L\ 
values; one can expect that this will influence studies in- 
volving rotation curve fitting and the TF-relation. 

In this paper I concentrate on the use of color pro- 
files as a diagnostic tool to investigate dust and stellar 
content of spiral galaxies. Other processes that may con- 
tribute to the broadband colors (e.g. emission from hot 
dust in the K passband) are ignored as they are expected 
to be small in most cases. The structure of this paper is as 
follows. In Sect. || the data set is described and the color 
profiles of the 86 spiral galaxies using the B, V, R, I and K 
passband data are presented. Section || describes the ex- 
tinction models and the stellar population models used in 
this paper and then compares these models to the data. In 
Sect. f|, I investigate the relation between the color prop- 
erties of the galaxies and the structural galaxy parameters 
derived in the previous papers of this series. Implications 
of the current measurements are discussed in Sect. [s] and 
the paper is summarized in Sect. 

2. The data 

In order to examine the parameters describing the global 
structure of spiral galaxies, 86 face-on systems were ob- 
served in the B, V, R, I, H and K passbands. A full de- 
scription of the observations and data reduction can be 



found in de Jong & van der Kruit (1994, hereafter Paper I) 
of which only the essentials are repeated here. The galaxies 
in this statistically complete sample of undisturbed spirals 



were selected from the UGC (Nilson 1973) to have red di- 
ameters of at least 2' and axis ratios larger than 0.625. 
The galaxies were imaged along the major axis with a 
GEC CCD on the lm Jacobus Kapteyn Telescope at La 
Palma in the B,V, R and I passbands and with a near- 
IR array on the United Kingdom Infra-Red Telescope at 
Hawaii in the H and K passbands. Standard reduction 
techniques were used to produce the images, which were 
calibrated using globular cluster standard star fields. The 
sky brightness was determined outside the galaxy in ar- 
eas free of stars and its uncertainty constitutes one of the 
main sources of error in the derived parameters. 

The ellipticity and position angle (PA) of each galaxy 
were determined from the R passband image at an outer 
isophote (typically at 24 i?-mag arcsec -2 ). The radial sur- 
face brightness profiles were determined for all passbands 
by calculating the average surface brightness on ellipti- 
cal annuli of increasing radius using the previously deter- 
mined ellipticity and PA. This method ensures that the 
luminosity profile reflects the average surface brightness 
at each radius, independent of passband. Free ellipse fit- 
ting at each surface brightness interval (e.g. Kent 1984 ) 
will give a disturbed representation of the radial luminos- 
ity distribution in face-on galaxies, due to spiral arms, 
bars and bright H n regions. The surface brightness pro- 



files were used to calculate the integrated luminosity of 
the galaxies. Internal and external comparisons showed 
that the derived parameters are well within the estimated 
errors. 

The decomposition of the light of the galaxies into its 
fundamental components (bulge, disk and sometimes a 
bar) is described by de Jong (1996a, hereafter Paper II). 
An exponential light distribution was assumed for both 
the bulge and the disk and these were fitted to the full 
2D image. An extensive error analysis of the determina- 
tion of the fundamental galaxy parameters was performed 
and this revealed that the dominant source of error is the 
uncertainty in the sky background. 

The color profiles of the galaxies were calculated by 
subtracting the radial surface brightness profiles of the dif- 
ferent passbands from one another. Some typical profiles 
are presented in Fig. [I , the profiles of all galaxies can be 
found in de Jong ( 1995 ) and are available in electronically 
readable format. The dashed lines indicate the maximum 
errors due to the uncertainty in the sky surface bright- 
ness. One should be cautious in interpreting the colors in 
the inner few seconds of arc, because the profiles were not 
corrected for the differences in seeing (Paper I) between 
the different passbands. 

A quick inspection shows that almost all galaxies show 
color gradients. They become bluer going radially out- 
ward, even when taking the sky background subtraction 
uncertainties into account. The color gradients extend over 
several disk scalelengths. Note that bulges leave no clear 
signature in the color profiles. From the color profiles alone 
one can not tell which part is bulge dominated and which 
part is disk dominated. 

The profiles in Fig. [l] are the observed profiles. The cor- 
rections needed to translate observed quantit ies into more 
physical quantities are discussed by de Jong ( 1996b , here- 
after Paper III). In the remainder of this paper, only the 
photometric observations are used, corrected for G alactic 
extinction using the precepts of Burstein and Heiles ( 1984 ) 
and the extinction curve of Rieke and Lebofsky ( 1985 ) . 



3. Color gradients 

The color gradients of Fig. [I] have been put on a common 
scale in Fig. ||, where I have plotted for all galaxies the 
average B-K color at each radius as function of the az- 
imuthally averaged R passband surface brightness at the 
corresponding radii. The galaxies are divided in four bins 
based on their morp hological type, using the RC3 (de Vau- 
couleurs et al. |l99l| ) type indices T (see also Papers I and 
III). There is a clear correlation between average surface 
brightness at a radius and the average color at that radius; 
the lower surface brightness regions are bluer. This indi- 
cates the relation between Hubble type, surface brightness 
and integrated color: since late-type galaxies have on av- 
erage a lower central surface brightness (Paper III) , they 
arc bluer. This is not the whole story, since for each mor- 
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Fig. 1. Observed color profiles for 8 out of the 86 galaxies. The dashed lines indicate the maximum error due wrong sky 
background subtraction, the arrow indicates one disk scalelength in the K passband. 
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Fig. 2. The average B-K colors of all galaxies at different radii as function of the azimuthally averaged R passband surface 
brightness measured at these radii. The galaxies are divided into the four indicated RC3 morphological type index (T) bins. 
The dashed lines give a common reference in all four bins, but have no physical meaning. The lines have a B-K color gradient 
of 1/7 mag per 7?- mag arcsec~ 2 . 



phological type at each surface brightness there is con- 
siderable scatter. Furthermore, even at the same average 
surface brightness, late-type galaxies are on the average 
bluer than early-type galaxies. 

The two most straightforward explanations for the 
color gradients are 1) radial changes in stellar popula- 
tions and 2) radial variations in reddening due to dust 
extinction. For both possibilities, I investigate a range of 
models to limit the acceptable parameter space. The ex- 
tinction models have a range in relative distributions of 
dust and stars. The colors of the stellar synthesis popula- 
tion models depend on the star formation history (SFH) 
and the metallicity of the stars. 

The colors and color gradients of the galaxies formed 
from the different passband combinations are correlated 
and the models should be fitted in a six-dimensional "pass- 
band space" . Predicting the right color gradient in one 
combination of passbands, but a wrong one in an other 
combination makes a model at best incomplete and there- 
fore undesirable. Color-color plots will be used to show as 



much information as possible in one plot. The H passband 
data are not shown, as the differences between H and K 
predicted by the models (both the population and the ex- 
tinction models) are smaller than the measurement errors. 
In the remainder of this section I first discuss the extinc- 
tion models and the stellar population synthesis models 
used in this paper and then compare the models with the 
data. 

3.1. Extinction models 

In this section, I present my new dust models, show the 
predicted luminosity profiles, color profiles and color-color 
diagrams, and compare the results with existing dust mod- 
els. 

3.1.1. Modeling dust effects 

Numerous researches have investigated the effects of dust 
extinction on the observed light distributions of galaxies. 
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The primary goal of most of the studies is to investigate 
the inclination depende nt eff ects of the total magnitude of 
galaxies (e.g. Huizinga 1994| ). In some studies the extinc- 
tion effects on the observed (exponential) light profile of 
galaxies is studied. The most detaile d are the Triplex mod- 
els by Disney, D avies fc Philli pps (1989 , hereafter DDP; 
see also Huizinga 1994 ; Evans |l994 ) . 

The effects of reddening on the observed colors and 
color profiles has been examined in a number of studies. 
The simplest model to predict the reddening of a galaxy 
is to use directly the standard (Galactic) extinction law, 
but this is of course a gross oversimplification. DDP have 
named this the Screen model, which has all dust placed 
between us and the galaxy. In reality the dust is mixed 
between the stars, so that on the near side of the galaxy 
a considerable fraction of stars will be only slightly ob- 
scured. For the same amount of dust, the observed red- 
dening is considerably less than predicted by the Screen 
model, especially since the most reddened stars are also 
the most obscured stars and therefore the ones that con- 
tribute less to the overall color of the system. 

As soon as the dust is mixed with the stars one has 
to take both absorption and scattering into account. In- 
tuitively one expects that for face-on galaxies at least as 
much light gets scattered into the line of sight as out of 
it, especially since there are more photons traveling in the 
plane of a galaxy which can be scattered into face-on di- 
rections than the other way around. As only the absorbed 
photons really disappear, it is better to use relative ab- 
sorption rather than relative extinction between different 
passbands to estimate reddening effects in face-on gal- 
axies. It is essential to incorporate both absorption and 
scattering into extinction models to make accurate pre- 
dictions of the effects of dust on colors and color gradients 
of galaxies. 

A number of studies have investigated the effect of 
red denin g on integ rated colors of galaxies (Bruzual et 
al. 1988; Witt et al. 1992 and references therein). In these 



studies scattering is included and stellar and dust dis- 
tributions are used that allow approximations to reduce 
computing time; e.g. Bruzual et al. use plane parallel dis- 
tributions and Witt et al. use spherically symmetric distri- 
butions. Color profiles produced by dust models are not 
often presented. Evans ( 1994 ) investigates the effects of 
extinction as function of radius in face-on galaxies for a 
non-scattering medium. Byun et al. (1994) also investigate 
the effects of dust on luminosity an d colo r profiles, using 
the method of Kylafis & Bahcall ( |1987| ). Their method 
includes first order scattering and approximates multiple 
scattering. The results of Byun et al. 1994 are compared 
with the results presented here in Sect. 3.1.2. 



To estimate to what extent the color gradients can be 
attributed to reddening by dust extinction, Monte Carlo 
simulations were made of light rays traveling through a 
dusty medium. The models are described in full detail in 
Appendix 



The distributions of stellar light and dust in these mod- 
els were described by exponential laws in both the radial 
and vertical directions. In the radial direction these dis- 
tributions were parameterised by the scalelength of the 
stars (/i s ) and the dust (/id), and in vertical direction by 
the scaleheight of stars (z s ) and dust (zd). In all models 
I used h s /z s — 10 and for simplicity no bulge component 
was added to the stellar light distribution. 

Since the effects of dust on the color profiles is the main 
interest of this study, absolute calibration of the amount 
of starlight is arbitrary. Only the relative effect of dust 
from one passband to the other is important. The amount 
of dust in the models is parameterised by the optical depth 
of a system, defined as the optical thickness due to 

dust absorption and scattering in the V passband through 
the disk from one pole to the other along the symmetry 
axis (Eq. (|At6|)). 
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Fig. 3. The relative extinction (ta/tv) and the relative ex- 
tinction ((1 — (i\)t\/tv) used in the dust models. 



Three dust properties were incorporated into the dust 
model to describe the wavelength dependent effects of 
dust extinction: the relative extinction (ta/tv), the albedo 
(ax), and the scattering asymmetry parameter (g\). The 
relati ve ex tinction was adopted from Rieke and Lebof- 
sky ( 1985 ); the othe r two parameters were drawn from 
Bruzual et al. ( 1988j ). These values are listed in Table [j] 
and Fig. || shows the relative extinction and the relative 
absorption ((1 — a\)r\/Tv). Note that for the optical pass- 
bands the change in relative absorbtion is much smaller 
than the change in relative extinction. 

Before presenting the model results, a few words of 
caution are in order. First, extragalactic dust properties 
are poorly known. There are only a few measurements 
of extinction laws in extragalactic systems (e.g. Knapen 
et al. 1991 ; Jansen et al. 1994 ) other than for the Magel- 
lanic Clouds (see Mathis 199C for references). All measure- 
ments seem to be consistent with the Galactic extinction 
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passband 


t\/t v 


a\ 


9x 


U 


1.531 


0.68 


0.67 


B 


1.324 


0.66 


0.59 


V 


1.000 


0.60 


0.50 


R 


0.748 


0.53 


0.40 


I 


0.482 


0.45 


0.29 


H 


0.175 


0.28 


0.04 


K 


0.112 


0.20 


0.00 



Table 1. Values used in the scattering model for the dust prop- 
erties relative extinction (t\/tv), albedo (a\) and scattering 
asymmetry (g\) as function of photometric passband. 



law, except for a few measurements in the Small Magel- 
lanic Cloud. It is well known that the Galactic extinction 
curve is not the same in all directions, but the one adopted 
here is appropriate for the diffuse interstellar medium (for 
discussion see Mathis 1990 ). The parameters a\ and g\ 
have never been measured in extragalactic systems and 
are poorly known even for our own Galaxy. The adopted 
values for these parameters stem, especially for the longer 
wavelengths, from model calculations. Still, no large varia- 
tions are expected in the extinction properties, unless the 
dust in other galaxies is made of totally different material 
(see also the discussion in Bruzual et al. 198S| ). 

As a second word of caution, the models presented 
here describe only smooth diffuse dust. The effects of non- 
homogeneous dust distributions should be considered. A 
large ensemble of optically thick clouds has only a red- 
dening effect if the clouds have a large filling factor, but 
such a configuration becomes comparable to the presented 
models with high r y. The reddening effect of a clumpy 
medium will be smaller than the effect predicted by the 
diffuse dust models for the same amount of dust, but the 
direction of the reddening vectors will be the same as long 
as the dust properties in the clouds are more or less the 
same as used here. If clouds are optically thick at all wave- 
lengths one has the case of gray dust and no color gradients 
at all. Model calculations using a clumpy dust medium in 



the absence of scattering are presented in Huizinga (1994, 
Chapter 5). 

As a final word of caution, a young population of stars 
probably has a smaller scaleheight than an old population 
of stars. It might be more appropriate to use a smaller 
stellar scaleheight in the blue than in the near-IR. The 
relative contributions from young and old populations are 
difficult to estimate however, and for simplicity one stel- 
lar scaleheight is used for all passbands. These models do 
not include the dust shells around the extremely luminous 
stars in the final stages of their life. Even though such 
shells will make these stars redder, they will not produce 
a radial effect (unless the shell properties depend on galac- 
tic radius). Effectively these shells will only make the total 
underlying population redder at all radii and they are of 
no further concern here. 



3.1.2. Resulting profiles 

Figure ^ shows luminosity and color profiles resulting from 
the Monte Carlo simulations. The luminosity profiles for 
the different passbands have been given an arbitrary off- 
set and the dust free cases of the B and the K passbands 
are indicated by the dashed lines. The color profiles have 
been plotted under the arbitrary assumption that the un- 
derlying stellar populations have color indices of zero in 
all passband combinations. The noise in the color profiles 
is due to the statistical processes inherent to Monte Carlo 
simulations. 

The luminosity profiles of the hd/h s = l models show 
deviations from the unobscured profiles only at the inner 
two scalelengths. These deviations are quite small except 
for the highest t .v values. The differences between the 
different z&j z & models are also quite small and are only 
apparent for the high Toy values. The gradients can be 
large in B-H and B-K, but are in general smaller than 
0.3 mag in the other color combinations. The color gra- 
dients are small in the wavelength range from the U to 
the R passband, because the absorption properties do not 
differ very much among these passbands. The change in 
scattering properties causes the differences in extinction 
in this wavelength range. 

The luminosity profiles for the hd/h s = 2-3 models are 
affected over several scalelengths by dust extinction. In 
fact, the hd/h s = 3, tq,v = 10-20 models are optically thick 
over almost the entire disk. The result is that the pro- 
files stay exponential, but with a lower surface brightness 
and a slightly different scalelength from the unobscured 
case. The color profiles show no color gradients for these 
models, only color offsets. The typical surface brightness 
is produced at t\ — 1 over the entire disk, and the color 
offsets reflect that different wavelengths probe different 
depths into the galaxy. 



are 



To my knowledge, the models of Byun et al. 1994 
the only models in the literature that have exponential 
light and dust distributions, and include scattering to cal- 
culate luminosity and color profiles. These models can be 
compared to the models presented, but only indirect, be- 
cause Byun et al. defined the optical depth of a system dif- 
ferently. They parameterize the optical depth of a system 
as the absorption in the V passband through the whole 
disk of a face-on galaxy along the s ymmetry axis, while 
here the extinction is used (Eq.( |A16| )). Furthermore, they 
use the Galactic extinction law to translate the absorption 
coefficient from one passband to another. Using Table [j] 
one can calculate that their ry(0) models corresponds to 
m y Toy~Tv(Q)/(l — ax) models, which is 2.9ry(0) for the 
B passband and 1.8ry(0) for the / passband. It is proba- 
bly most meaningful to compare the tq y = 20, z&j z s — 0.3, 
hd/h s = 1 B passband profile of Fig. |j with the bulgeless 
(BT0.0), face-on r(0) = 5.0 profile of their Fig. 7. The cen- 
tral extinction of slightly more than 1 mag and the gen- 
eral shape of the luminosity profile (which is unaffected by 
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extinction for radii larger than 2-3 scalelengths) are com- 
parable. Their B—I color profiles are physically not very 
plausible, because they use the Galactic extinction curve 
instead of an absorption curve to translate their absorp- 
tion coefficients from one passband to another. A compar- 
ison of the color profiles is therefore not meaningful. Their 
luminosity profiles can be used, but note that their Ty(0) 
should be divided by (l-ay) to get the extinction suffered 
by a point source behind the galaxy, as done here. 

What is the range of plausible model parameter val- 
ues? Quite high Toy and/or hd/h s values are needed to 
explain the observed B-K color gradients of one mag- 
nitude over five scalelengths (Fig. [j]) by dust reddening 
alone. On the other hand, models with equal scaleheight 
for dust and stars need an additional dust component, be- 
cause these models do not produce a clear dust lane in 
edge-on galaxies. It is also unlikely on dynamical grounds 
that the dissipational dust and the dissipationless stars 
have the same scaleheight. Kylafis & Bahcall ( 1987 ) find 
a dust-to-star scaleheight ratio of 0.4 in their best fitting 
model of edge-on galaxy NGC 891 and the dominant dust 
component is expected to have a Zd/z s ratio between 0.3 
and 0.5. 

The high h^/h s models are favored by Valentijn ( 1990| , 



1994), who concluded from inclination tests that Sb-Sc 
galaxies have a Tg ~ 1 through the disk at Z?25- This 
extinction at about 3-4 stellar scalelengths translates to 
Toy ~ 20 models, if the dust density is distributed expo- 
nentially with hd/h s — l. The edge-on extinction from the 
center out to 3-4 stellar scalelengths gives at least ry « 50 
in such a model, which is in conflict with observations of 
edge-on galaxies. A few edge-on galaxies have been im- 
ag ed in the near-IR indicating Ay ~ 8-10 (Wainscoat et 
al. 1989 ; Aoki et al. 1991 ), and the Galacti c cent er can be 
seen in the K passband (Rieke & Lebofsky 1985] , Ay ~ 30 
and Ak ~ 3). Thus Toy — 20 models represent extremely 
dusty galaxies, and certainly no galaxies with To ; y>20 are 
expected. If one increases hd/h s to 3, Toy has to be 2-4 
to get Ty = 1 through the disk at 3 stellar scalelengths. 
This then gives an edge-on Ty = 20-40 from the center out 
to 3 stellar scalelengths when using large Zd/z s values, 
which are the most favorable for these models. Therefore, 
the edge-on extinction values of the Toy = 2-4, hd/h s = 3 
models are marginally consistent with the observations, 
but these models do not produce very large color gradi- 
ents. 



indicated Toy values were multiplied by (tv/t\)(1 — a\)/2 
(Table 0). This is equivalent to using an absorption law 
rather than an extinction law between the different pass- 
bands. The factor 2 arises because the DDP models are 
characterized by the optical depth from the galaxy cen- 
ter to the pole and not by the optical depth through the 
whole disk. 

Comparing the results from the Monte Carlo simula- 
tions with the DDP models, one can see that the models 
agree remarkably well for high optical depths. The intu- 
itive idea that just as many photons are scattered out 
of the line of sight as are scattered in seems correct for 
face-on galaxies. Photons are only lost due to absorption. 
For low optical depths the reddening almost completely 
disappears. Once a photon gets scattered into the line of 
sight, the chances of it getting absorbed or scattered again 
are minimal; even bluing instead of reddening can occur. 
Even though the amount of reddening is a strong function 
of the dust configuration, it seems that the direction of the 
reddening vector is largely dependent on the dust proper- 
ties. Obviously the reddening produced by these models is 
different from the reddening produced by a Screen model 
with the Galactic extinction law (also indicated in Fig. ^J) . 

In conclusion, dust can produce color gradients in face- 
on galaxies, but this requires quite high central optical 
depths and preferably long dust scalelengths. The redden- 
ing vectors of realistic dust models that include both ab- 
sorption and scattering arc completely different from the 
often-used Screen model extinction models. 

3.2. Evolutionary stellar population synthesis models 

Ever since the invention o f the concept of stellar popu- 
lations in galaxies (Baade |l944 ), numerous models have 
been made to predict the integrated light properties of 
such populations and thus of galaxies as a whole. Ini- 
tially the empirical approach was often followed, in which 
the different contributions of the stellar populations are 
added to match the observed galaxy SED. Later, knowl- 
edge about initial conditions and stellar evolution were 
added to make evolutionary synthesis models. This sec- 
tion contains a brief description of population synthesis 
models, concentrating on the synthesis models used here, 
followed by a description of the effects that star formation 
history (SFH) and metallicity have on the colors produced 
by synthesis models. 



3.1.3. Resulting color-color diagrams 



3.2.1. Modeling stellar populations 



Four color-color plots of dust models are presented in 
Fig. ^ Because extinction is a relative measurement, the 
zero-point can be chosen freely in these plots; only the 
shapes of the profiles are fixed. The solid lines show the 
results from the Monte Carlo simulations. The dotted lines 
are the result of the Triplex models of DDP. To calculate 
the To for the DDP models in the different passbands, the 



In recent years a large number of synthesis methods have 
appeared in the literature, all of which require many in- 
put parameters (Tinsley 1980| ; Renzini & Buzzoni 1986; 



Worthey 1994 and references therein). The simplest mod- 
els use a initial mass function (IMF) to create a single 
burst of stars, whose evolution in time is then followed. 
Slightly more complicated models describe the star for- 
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mation rate (SFR) in time. In the most complicated mod- 
els, the gas, stellar, and chemical evolution are linked and 
described in a self-consistent way. This last type of model 
has not been used here, since to date only a small range of 
SFHs have been investigated with these models. The de- 
scription of the evolution alone does not produce an SED 
and therefore the models are linked to a stellar library to 



calculate the evolution in time of the integrated passband 
fluxes or of the integrated spectrum. 

The results from the models in the literature are not 
all in agreement. The disagreements arise mainly from dif- 
ferences in the treatment of the late stages of stellar evo- 
lution. The models do agree on the two main parameters 
determining the integrated colors of a synthesized galaxy. 
First of all the colors are strongly determined by the col- 
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ors of the youngest population, thus by the SFH, and sec- 
ondly the colors are considerably affected by the metallic- 
ities of the populations. In Sect. |l] it was noted that both 
these SFH and metallicity changes have been observed 
on a radial scale in spiral galaxies. Furthermore, the ra- 
dial age and metallicity gradients in our own Galaxy are 
well known and have been exten sively studied (Gilmore 
et al. 198£ ; Matteucci 1989 , 1992 and references therein). 
Synthesis models incorporating both age and metallicity 
effects are needed in the comparisons with spiral galaxy 
data. 

The population synthesis models of Bruzual & Char- 
lot Q1996| , BC96 models hereafter, see also Chariot 



& Bruzu al [1991 Bruzual & Chariot ( |1993| ) and of 
Worthey ( 1994 , W94 models hereafter) are used in the 
remainder of this paper; they are shown in the color- 
color diagram of Fig. g. The BC96 models used here are 
based on the is ochrone tracks of the Padova group (Bres- 
san et al. 1993) and on an empirical stellar flux library. 
The BC96 isochrone synthesis approach makes calcula- 
tion of the very early stages of evolution of a popula- 
tion possible. The W94 mo dels a re constructed from the 
isochrones of VandenBerg (1985) and the Revised Yale 
Isochrones (Green et al. 1987) and use a theoretical stel- 
lar flux library. A different approach was followed in the 
BC96 and W94 models to calculate the integrated spectra 
of an evolving stellar population, but the main difference 



of interest here is the regions of age-metallicty parameter 
space that were investigated. The BC96 models were only 
calculated for solar metallicity, but give colors of popula- 
tions as young as 1.26xl0 5 yr. The W94 models span a 
wide range in metallicity, but the youngest population is 
1.5 Gyr. 

The models used here were calculated with the stan- 
dard Salpeter IMF (Salpeter |l955| ), with lower mass cut- 
off at 0.1 Mq and upper mass cutoff at 2Mq (W94) or 
100 Mq (BC96). The use of other solar neighborhood type 
IMFs or other cutoffs has only small effects compared 
to the main factors determining the colors of a synthe- 
sized population, namely age (or SFH) and metallicity. 
The main result will be some shifts of the colors in color- 
color space. Using IMFs deviating strongly from the solar 
neighbourhood IMF will result in substantial changes, but 
without strong evidence these types of IMFs are hard to 
justify. 

A comparison study of several stellar population syn- 
thesis mod els ha s been performed by Chariot, Worthey 
& Bressan ( 1995 , hereafter CWB). They find substantial 
discrepancies among the models, especially in V-K, which 
are larger than the typical observational errors for gal- 
axies. They conclude that age determinations from galaxy 
colors are about 35% uncertain at a given metallicity and 
25% uncertain in metallicity at a given age, even for a 
single burst population. The main source of disagreement 
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stems from differences in the stellar evolution prescription, 
with a smaller component resulting from the spectral cal- 
ibrations. For comparison, to see the effect of other stel- 
lar evolutionary tracks, de Jong ( 1995 , 1996c ) presents a 
similar color-color analysis as presented here, but using 
Bruzual & Chariot (|l993| ) with Maeder & Meynet ( |l99l[ ) 
evolutionairy tracks instead of BC96 with Padova tracks. 

The W94 and B96 models were chosen to cover a large 
age-metallicity space, but given the uncertainties in the 
models, any other set of up-to-date set of models could 
have been chosen (e.g. Bres san, Chiosi & Fagotto 1993| ; 
von Alvensleben & Gerhard 1994 ). The particular choice 
of models will not change the main results of this study. 
The model uncertainties will mainly result in the uncer- 
tainty of absolute calibration of colors, but the relative 
color trends are reasonably correct. With the current state 
of modeling, I will only indicate which regions in gal- 
axies and what type of galaxies are younger/older or have 
higher/lower metallicity. 

3.2.2. Star formation history in color-color diagrams 

Three evolutionary tracks of the BC96 models arc shown 
in Fig. 1^. The simplest is the single burst model, in which 
the color evolution of one initial starburst is followed in 
time. In the other two tracks this single burst model has 
been convolved in time to yield color evolution for differ- 
ent SFHs. In one model an exponentially declining SFR 
with a time scale of 5 Gyr was used, in the other model 
the SFR was held constant. These models show the im- 
portance of the very early stages of stellar evolution. The 
constant SFR model at 17 Gyr is as blue in B-V and B- 
I as the exponentially declining model at 8 Gyr and the 
single burst model at 1.5 Gyr! Still, a solar metallicity star- 
burst cannot produce colors to the right of the indicated 
single burst line in Fig. |^, because the colors of very young 
populations (age < 1.5 Gyr, not shown) all lie to the left 
of (or in the direction of) the single burst trend. 

3.2.3. Age and metallicity in color-color diagrams 

Of the W94 models, only the single burst models of differ- 
ent metallicities are shown in Fig. ^. The synthesis method 
of W94 prevents calculation of the very early evolution 
stages of a stellar population. Young populations are espe- 
cially hard to synthesize for the lower metallicities, simply 
because there are no young, low metallicity stars in the so- 
lar neighborhood that can be used as input stars for the 
models. 

The lower metallicity populations of the W94 models 
are clearly bluer in all color combinations. Age and metal- 
licity are not complete degenerate. The offset in optical- 
near-IR colors is slightly larger than the offset in optical- 
optical colors. A low metallicity system can be recognized 
by its blue V— K or R-K color with respect to its B-V 
and B-I colors. When mixing populations of different age 



and metallicity the model grid points can be more or less 
added as vectors and then age and metallicity are degen- 
erate. An individual galaxy can never be pinpointed to a 
certain age and metallicity using broadband colors. 

The single burst BC96 and W94 models are clearly 
offset from each other. The solar metallicity BC96 model 
lies somewhere in between the [Fe/H] = —0.25 and [Fe/H] 
= 0.0 W94 models. These are the kinds of differences be- 
tween models that were studied in detail by CWB, re- 
sulting mainly from uncertainties in stellar evolutionairy 
tracks. The relative trends agree reasonably well, confirm- 
ing that it is allowed to look the relative trends in color 
space. The BC96 models will be mainly used to investi- 
gate trends in SFH, because contrary to the W94 models, 
they incorporate the very early stages of stellar evolution. 
To compare with the data, I will connect the 8 Gyr points 
of the different SFH BC96 models and do the same for the 
17 Gyr points. These connected data points will indicate 
the color trends for equally old populations with different 
SFHs. The W94 models have to be used to look at the ef- 
fects caused by metallicity, because different metallicities 
are not available in the BC96 models. 

3.3. Color gradients; measurements versus models 

In this section I first use color-color diagrams to display 
the observations. I then try to explain the observed color 
gradients by comparing these measurements with the dust 
models, the different SFH synthesis models, and finally the 
full population synthesis models, that incorporate both 
age and metallicity effects. 

3.3.1. The measurements in color-color diagrams 

The color-color profiles of the galaxies are presented in 
Figs |?] and ||. The data were smoothed to reduce the noise 
in the profiles. The first data point (at the open circle) 
is the average over the inner half scalelength (in the K 
passband) of the luminosity profiles, the other points are 
averages moving outward in steps of one K scalelength. 
The uncertainty of the inner point is dominated by the 
zero-point uncertainty of the calibration and is indicated 
at the bottom of the top-left panels of Figs and ||. In 
the blue direction of the color-color profiles, lower surface 
brightnesses are traced (see Fig. ||) and errors are domi- 
nated by sky background uncertainties. Strange kinks at 
the blue ends of the profiles should thus not be trusted. 

The profiles of galaxies with T<6 are confined to a 
small region in the color-color plots. The colors between 
different passband combinations are strongly correlated. 
The scatter is slightly larger than the average zero-point 
error. The central colors of the galaxies become on average 
a little bit bluer going from T=0 to T=6, but they follow 
the main trend. 

The galaxies with morphological classification T>6 
clearly deviate from the main trend. Their central col- 
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Fig. 7. The B—V versus R—K color-color plots of the program galaxies divided into the four morphological type bins indicated 
at the bottom-right of the panels. The centers of the galaxies are indicated by the open circles; the thin lines follow the color 
profiles in radial direction in steps of 1 K passband scalelength. At the b otto m-left of the 0<T<3 panel, the typical rms error in 
the zero-point calibration is indicated. Some of the dust models of Sect. 3.1 are plotted in the top-left of the panels. Th e th ick 
dashed lines in the center of the panels connect the 17Gyr and the 8Gyr points of the BC96 models presented in Sect. 3.2. In 



the 5<T<6 panel the 12Gyr single burst W94 models for different metallicities are connected, the marks indicate the [Fe/H] 
values. The thick solid lines in the 6<T<10 panels represent the W94 models. Points of equal metallicity, but different age are 
connected. 



ors range from the normal red to extreme blue, even bluer 
than the bluest outer parts of the earlier type galaxies. 
The spread in the color-color diagram is also significantly 
larger for the late-type galaxies. Some of the late-type gal- 
axies have very blue V—K and R—K colors for their B—V 



and B—I colors, especially when compared to the earlier 
types. 

When comparing models with the data one should re- 
alize that as soon as a model for a galaxy has been chosen, 
it should be applied to all color combinations. In particu- 
lar, the same model should be used in both Figs fj] and |[ 
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Fig. 8. Same as Fig. [?], but for B-I versus V-K. 



It is tempting to propose one single model for all galaxies, 
because the profiles are confined to a small region in the 
diagrams. We only need to explain the offsets from the 
main trends with additional parameters. 

3.3.2. Measurements versus dust models 

The reddening profiles produced by the dust models are 
indicated on the left in the panels of Figs fj] and ||. The 
color of the underlying stellar population is arbitrary and 
thus the dust profile can be placed anywhere in the di- 
agram. The dust models have a distinct direction in the 
color-color diagrams independent of the dust configura- 



tion, as explained in Sect. 3.1. This direction is clearly 
different from the general trend of the data and therefore 
the whole gradient cannot be produced by the dust red- 
dening alone. A small fraction of the color gradients could 
be due to dust reddening, but an additional component is 
needed to explain the full gradient. 

This does not mean that there could not be large 
amounts of dust, but rather that the color gradients are 
not mainly caused by dust reddening. If the dust is not 
diffuse, but strongly clumped into clouds the amount of 
reddening is strongly reduced. It could be that a large 
fraction of the most luminous stars is embedded in dust 
clouds. This will not induce a color gradient or an inch- 
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nation dependent extinction effect, but will give the total 
color profile an offset in the general direction of the calcu- 
late d dus t models. The "dusty nucleus" models of Witt et 
al. ( 1992 ) give an even better indication of the expected 
offset vector if the luminous stars are embedded in dust 
clouds. 

3.3.3. Measurements versus metallicity effects 

To what extent can metallicity differences in the stellar 
populations account for the color gradients? The 12-Gyr- 
old W94 single burst models for different metallicities are 
connected in the 5<T<6 panels of Figs ^| and |[ Models 
at other ages follow the same direction in color space. The 
single-age, different-metallicity model does not match the 
data for most of the galaxies. Again another component 
is needed to explain the color gradients; radial metallicity 
differences alone are not sufficient. It is important to note 
that population models cannot be arbitrarily shifted, they 
predict fixed colors for a given SFH and metallicity. 

3.3.4. Measurements versus star formation history 

The BC96 models are indicated by dashed lines in Figs [7] 
and ||. The 8 Gyr points of the models have been connected 
as well as the 17 Gyr points. The reddest ends of these 
lines indicate single burst models, the bluest ends repre- 
sent constant SFR models. In between is a model with an 
exponentially decreasing SFR with a time scale of 5 Gyr. 
The color vector of the T<6 systems seem to be reason- 
ably well matched by the BC96 models, but most galaxies 
are offset to the blue in the V—K and R-K colors. The 
offset is the smallest for the 8 Gyr model, but then most of 
the galaxy centers are even redder than predicted by the 
single burst model. Because we know that galaxies still 
have star formation in their centers, this is an unlikely sit- 
uation. As mentioned before, model uncertainties allow for 
some shifts in these diagrams, but the most likely conclu- 
sion has to be that although SFH variations as function of 
radius are a good driving force for radial color gradients, 
alone they cannot explain the full color gradients. This is 
especially true for the galaxies with T>6, which certainly 
have V—K and R-K colors too blue to be explained by 
the solar metallicity BC96 models. 

3.3.5. Measurements versus age and metallicity effects 

The W94 models in the 6<T<10 panel indicate that the 
very blue galaxies in this panel can be described very well 
by low-metallicity population synthesis models. As with 
the earlier type galaxies, the radial color trends are rea- 
sonably well described by age differences, but at all radii 
lower metallicities are needed for most of the galaxies. A 
number of them are so blue in all color combinations that 
their stellar components must be young and of low metal- 
licity. 



Galaxies are known to have radial metallicity gra- 
dients in th eir current gas content (Villa -Cost as & Ed- 
munds 1992 hereafter VE; Zaritsky et al. 1994 , hereafter 
ZKH), and have radial SFRs that are not linearly cor- 
related with their radial stellar surface brightness, which 
means they do n ot ha ve one SFH as function of radius 
(Ryder & Dopita 1994 ). As long as there are no consis- 
tent stellar population synthesis models that incorporate 
very young stellar evolutionary stages at all metallicities, 
it is difficult to make quantitative statements about the 
observed colors and color gradients of the galaxies. 

Limits can be set on the models using the metallic- 
ity measurements collected by VE and ZKH. Recall that 
metallicity measurements yield the current gas metallic- 
ities in Hn regions, so that the underlying stellar com- 
ponent could have completely different metallicity values. 
For galaxies T<6 the 12+log(0/H) values run from ~9.3 
in the center to ~8.6 at R25. The O/H indices of later 
type galaxies are a few tenths lower, from about 8.9 to 8.2, 
but with more or less the same gradient across the disk. 
Using log(O/H) Q ~-3.08 (Grevesse & Anders |98|) and 
[0/Fe]~ 0.1-0.5 (e.g. Wyse & Gilmore |l988|) the O/H 
values can be transformed to [Fe/H] values and used with 
the W94 models. The [Fe/H] values run from about a cen- 
tral 0.2 to -0.6 in the outer regions in galaxies with T<6, 
and from approximately -0.1 to -1.1 in the later type gal- 
axies. 

The metallicities just calculated should not be taken 
too literally in the comparison of the W94 models with the 
data in Figs [7] and ||. As noted before, models should be 
used to indicate trends in color-color space but cannot be 
expected to give absolute colors in an individual galaxy. 
Still, the B-V and B—I colors are far too red at each ra- 
dius using the 12 Gyr W94 models in the metallicity range 
determined from the Hil regions ([Fc/H]=0.2 to -0.6). 
This must mean that either the underlying stellar popula- 
tions have a much lower average metallicity than the sur- 
rounding gas or that at each radius a much younger stellar 
population is present, making the average age lower than 
12 Gyr. In the center, the contribution of the young stars 
cannot be very large, because the color vector of SFH (in- 
dicated for solar metallicity by the BC96 models) is in the 
wrong direction from the (t=12Gyr, [Fe/H] =0.25) point 
in the 5<T<6 panel. Therefore the center probably con- 
tains a mix of populations of different metallicities, with 
an average metallicity far lower than the current metallic- 
ity of the surrounding gas ([Fe/H] ~0.2), but on average 
higher than the gas metallicity of the outer regions. To- 
wards the outer parts of galaxies, the SFH color vector 
gives an excellent description of the observed color trend. 
In the outer regions it is unlikely that the metallicity of 
the stars is much higher than that of the gas for all gal- 
axies, and the metallicity of the stars must be of the order 
[Fe/H] =-0.6 (-1.1 for the late types) or less. Thus the 
color gradients are probably driven by a combination of 
SFH and metallicity differences as function of radius. The 



R.S. de Jong: Near-IR and optical broadband observations of 86 spirals. IV Color profiles 



17 



outer parts of galaxies are clearly younger on average than 
the central regions. If part of the color gradients is also 
caused by reddening, the age effect must even be larger, 
because this is the only effect that has a color vector that 
can compensate the dust color vector. 

In conclusion, the centers of spiral galaxies with T<6 
contain probably a relatively old population of stars, with 
a range in metallicities. The populations in the outer re- 
gions are on average much younger and their metallicity 
is probably lower, as the gas metallicity is much lower 
than the metallicity needed to explain the central col- 
ors. Overall, later type galaxies have lower metallicity, 
and a number of them are dominated by very young, 
low- metallicity populations. The observed color gradients 
cannot be caused by reddening alone, if the dust proper- 
ties used in the models are more or less correct. Likewise, 
metallicity gradients cannot cause the color gradients on 
their own, the color vectors are in the same direction as 
the dust model vectors, but as metallicity gradients are 
observed in the gas, it is likely they are present in the 
stars as well. Age gradients across the disk have to be 
even larger if reddening or metallicity gradients are im- 
portant. It should be noted that the population models 
can predict both the right colors and the right color gradi- 
ents within one system of models, while the dust models 
can at best explain only the color gradients. 

4. Colors and the structural galaxy parameters 

Most of the information that can be extracted from the 
galaxy colors of this data set is contained in the previ- 
ous sections, but for most data sets such detailed radial 
color information is not available. In order to allow com- 
parisons with these other data sets, a number of relation- 
ships between colors and fundamental galaxy parameters 
are shown in this section. 

In the previous section I argued that color gradients de- 
rived from different passband combinations are correlated 
in such a way that stellar population differences seem to 
be the most reasonable explanation of the phenomenon. 
Figure || suggests that the slope of the color gradient has 
a universal value for all galaxies, but this is an oversim- 
plification. 

The change in scalelength (h x ) as function of passband 
can be used to parameterise color gradients in the disk. 
The bulge/disk decomposition technique and the method 
used to determine the central surface brightnesses and 
scalelengths for the current data set were described in 
Paper II. A trivial recalculation shows that each axis of 
Fig. ^ indicates approximately the color change per scale- 
length in the relevant passband combinations. This figure 
illustrates again that color gradients are correlated, but 
also that there is no universal value for the gradient for 
all galaxies. The values range from about 0.1 to -0.8 mag 
per scalelength in B-K. 
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Fig. 9. The difference in scalelength between the different pass- 
bands. Only points with errors smaller than 0.15 are plotted. 
Different symbols are used to denote the indicated morpholog- 
ical types. 

The change in scalelength was used to find correla- 
tions between the steepness of the color gradient and 
other structural galaxy parameters. A large number of 
structural parameters were investigated, but none of them 
showed a correlation with color gradient. The parame- 
ters investigated include: inclination, morphological type, 
central surface brightness, scalelength, bulge-to-disk ra- 
tio, integrated magnitude, integrated colors, bulge color, 
bar versus non-bar, H I and CO fluxes, far-infrared IRAS 
fluxes and colors, group membership and rotation velocity. 
Fluxes at the different wavelengths normalized by area or 
integrated K passband flux were also used in the correla- 
tions, but with a negative result. A weak correlation was 
found only between the steepness of the color gradient and 
the central color of the disk. 

In the literature integrated colors of galaxies are usu- 
ally used to determine galaxy properties. The integrated 
colors as function of type are presented in Table || for this 
galaxy sample and the B—I colors are plotted in Fig. [l0| 
There is a clear correlation between type and color, but 
the scatter is large. The integrated color of a galaxy is 
dominated by the color of the central region and even 
though color correlates with (central) surface brightnesses 
(Fig. ||), one should note that each morphological type 
comes in a range of central surface brightnesses (Paper III) 
which then explains the large scatter in the integrated col- 
ors. It is better to determine and compare the colors of 
galaxies at a fixed isophote when looking for correlations. 

In Fig. [ll], the central bulge and disk colors of the 
galaxies are compared. The colors of bulge and disk are 
clearly correlated. This could be expected, as most color 
profiles show no clear changes in color gradient in the 
bulge region. Excluding the three deviant points with 
bulge B-K colors bluer than 2 mag, the bulge is on av- 
erage 0.14±0.55mag redder in B-K than the central disk 
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Fig. 10. The Galactic reddening corrected integrated B—I col- 
ors of the galaxies as function of morphological type index. The 
crosses show the values averaged over the bins indicated by the 
horizontal bars. The vertical bars are the standard deviations 
on the mean values. Only the galaxies with an error of less 
than 0.5 mag in their color were used. 
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Fig. 11. The central surface brightness B-K color of the disk 
versus the effective surface brightness B-K color of the bulge. 
Different symbols are used to denote the indicated morpholog- 
ical type ranges. The dashed line indicates the line of equality. 



color, which means that the stellar populations probably 
do not differ very much. 

5. Discussion 

An important consequence of the color differences in and 
among galaxies is the implied change in the M/L\ val- 
ues. The comparison of the data with the dust and stellar 
population models indicates that the color gradients in 
the galaxies of this sample result mainly from population 
changes and that dust reddening only plays a minor role. 
The M/L\ values corresponding to the synthesis models 
presented in Figs [l] and || are listed in Table ||. 



In every passband, young populations have much lower 
M/L\ values than old populations. Young populations 
still contain very massive stars, which are very luminous 
for their mass, but expend their energy quickly. Continu- 
ally adding in young populations, as done in the exponen- 
tially declining and constant SFR BC96 models, decreases 
the M/L\ values drastically. Young massive stars are blue 
so that this effect is most pronounced in the B passband. 
The change in M/Lb is about a factor of 6-11 between the 
2 and 17Gyr for both the BC96 and the W94 models and 
M/Lb changes a factor of 3-5 between the single burst 
and the constant SFR BC96 models. 

Metallicity has an entirely different effect on the M/L\ 
ratios. The M/Lb values increase with metallicity, while 
the M/Lk values decrease with metallicity. The turnover 
point is somewhere between the / and the J passband. 
This was the motivation for W94 to recommend the I 
passband for standard candle work and for studies of M / L 
in galaxies. 

The recommendation of the I passband is not entirely 
obvious, because the choice of optimum passband depends 
on whether one expects extinction, age or metallicity to 
have the largest influence on the photometry. If extinction 
is expected to play a role, the K passband should be used, 
as the extinction in the K passband is ^4 times less than 
the extinction in the I passband. The K passband should 
also be preferred if differences in SFR are expected to be 
important; see for instance the change in M/L\ ratios of 
the different BC96 models. One should turn to the I and 
J passbands only when the metallicity differences among 
the different objects are large. 

What influence do the current observations have on 
studies depending on M / L\ ratios? I address two issues 
here: rotation curve fitting and the TF-relation. 

The principle of rotation curve fitting is simple. One 
tries to explain the distribution of the dynamical mass of 
a galaxy by assigning masses to its known "luminous" 
(at whatever wavelength) components. The dynamical 
mass distributions can be determined from optical rota- 



tion curves along the m ajor axis (e.g. Rubin et al. [1985 



Mathewson et al. 1992 ) or, more sophisticatedly, fro m H I 
or other t wo-di mension al vel ocity fields (e.g. Bosma 1978; 



Begeman 1987 ; Broeils 1992 ). The mass assignment to the 
Hi and other gas components is relatively straightfor- 
ward, but the mass assignment to the stellar components 
has proven troublesome. Generally, M/L values have been 
assigned to the different stellar components (bulge and 
disk) using the maximum disk hypothesis (van Albada 
et al. 1985; van Albada & Sancisi 1986). This exercise re- 
vealed the "missing light" problem: the dynamical mass 
of a galaxy is much larger than the maximum "luminous" 
mass, with the main discrepancy in the outer regions. 

In the game of rotation curve fitting, one often encoun- 
ters the use of B or at best R passband luminosity pro- 
files. Table |^ clearly illustrates the danger of using these 
passbands, especially if one considers the color gradients 
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RC3 type 


# 


B-V 


# 


B-R 


# 


B-I 


# 


B-H 


# 


B-K 


0<T<2 


7 


0.81 ± 0.19 


10 


1.35 ± 0.22 


8 


1.83 ± 0.31 


6 


3.46 ± 0.21 


10 


3.75 ± 0.28 


2<T<4 


14 


0.74 ± 0.13 


26 


1.20 ± 0.19 


24 


1.78 ± 0.23 


6 


3.14 ± 0.34 


26 


3.53 ± 0.28 


4<T<6 


19 


0.67 ± 0.15 


19 


1.20 ± 0.13 


20 


1.75 ± 0.10 


10 


3.28 ± 0.14 


21 


3.51 ± 0.26 


6<T<8 


11 


0.59 ± 0.21 


11 


1.06 ± 0.26 


12 


1.61 ± 0.23 


3 


2.96 ± 0.62 


8 


3.18 ± 0.45 


8<T<10 


4 


0.69 ± 0.14 


5 


1.10 ± 0.10 


5 


1.44 ± 0.24 







4 


2.99 ± 0.09 



Table 2. The integrated colors for different morphological type ranges using only the photometric observations of Paper I with 
errors in the colors less than 0.5 mag. The # lists the number of galaxies in each subsample. 



Gyr 



BC models 
s.b. exp. cnst. 



W94 
-2 -1 



models, [Fe/H] = 
-0.5 -0.25 0.0 



0.25 



2 0.90 

5 3.13 

8 4.33 

12 5.70 

17 9.91 



0.82 
1.02 
1.42 
2.25 
4.02 



B 
0.30 
0.67 
1.02 
1.43 
1.88 



passband 



2.60 3.39 
3.71 4.81 
5.01 6.45 



1.52 1.73 

- 3.32 3.95 

4.16 5.19 5.96 

6.26 7.41 8.70 

8.53 12.28 10.54 



2.42 
5.26 
8.05 
11.95 
17.15 



I passband 



2 1.27 

5 2.26 

8 2.82 

12 3.23 

17 4.71 



1.25 
1.29 
1.57 
2.07 
2.95 



0.79 
0.91 
1.28 
1.68 
2.13 



2.34 2.53 
3.04 3.30 
3.84 4.10 



2.45 
3.50 
4.39 



1.09 
1.86 
2.67 
3.58 
4.73 



1.03 
2.00 
2.77 
3.60 
4.74 



1.24 
2.21 
3.08 
4.17 
5.54 



2 0.43 

5 0.80 

8 1.13 

12 1.15 

17 1.54 



0.65 
0.62 
0.74 
0.87 
1.13 



K 
0.24 
0.44 
0.62 
0.78 
0.93 



passband 



1.70 1.50 
2.10 1.88 
2.59 2.21 



- 0.52 0.42 

- 0.78 0.72 
1.16 1.12 0.86 
1.81 1.35 1.06 
1.90 1.69 1.29 



0.42 
0.60 
0.75 
0.87 
1.04 



Table 3. The M /L\ values in solar units for the models pre- 
sented in Section 3.2. For the BC models single burst models 

and constant 



(s.b.), exponential declining SFH models (exp. 
SFR models (cnst.) are listed. The single burst models of W94 
are listed for several different metallicities. 



observed here. Independent of whether the observed color 
gradients are caused by age or metallicity gradients (or 
by reddening for that matter), the M/Lb ratios will be 
much higher in the center than in the outer regions. Con- 
sequently the "missing light" discrepancy between inner 
and outer regions is even larger than estimated by the use 
of the B passband profiles together with a constant M/L. 
So what is the optimum passband to be used for rotation 
curve fitting? 



Following the conclusions of Sect. p.3| , let us assume 
that the central region of a T<6 galaxy consists of old 
stellar populations with a range in relatively high metal- 
licities, say on average t = 12 Gyr and [Fe/H] = 0. It is not 
likely that the metallicity of the stars in the outer re- 
gions is higher than that of the gas. The outer popula- 
tions are younger than the inner populations and let us as- 
sume that their average parameters are [Fe/H] =-0.5 and 
t = 8Gyr. Table || shows that M/Lb changes by a factor 
2.09 for those two populations, M/Lj by a factor 1.47, and 



M I Lk by a factor 0.91. Repeating this exercise for lower- 
metallicity late-type systems gives similar results. As long 
as the outer regions of galaxies are younger than the in- 
ner regions (and the example used was not very extreme) , 
the K passband is the optimum choice for rotation curve 
fitting. This is especially true if extinction plays a role, 
which is expected to be the case for the highly inclined 
galaxies normally used for rotation curve fitting. 

The Mj L effects on integrated magnitudes (as used in 
the TF-relation) are less trivial. In Sect. |l| it was argued 
that the integrated colors are dominated by the central 
colors of galaxies. In Figs, ^and |§|the open circles indicate 
the central colors of the galaxies and are therefore repre- 
sentative of the integrated colors. For the T<6 galaxies, 
the central colors follow the same trend as the color gra- 
dients within galaxies (as might be expected from Fig. ||) 
and therefore the same argument as for rotation curve 
fitting can be applied to recommend the K passband for 
TF-relation work. Note that the differences in central color 
in the bins of T<6 are much smaller than the differences 
within the galaxies themselves and a small color correction 
term should be sufficient to translate all galaxies to a com- 
mon M/L scale. The central colors of T>6 galaxies show 
a different distribution in color-color space. They come in 
a wide range of ages and metallicities and their M/Lk 
values can easily differ by a factor of two and by a factor 
of four in M/Lj. This would introduce an uncertainty of 
~0.75 if-mag or ~1.5 /-mag in the TF-relation respec- 
tively, if one would simply assume that the TF-relation 
is tracing the connection between luminous mass and dy- 
namical stellar mass. A two-color correction might then 
be needed to reduce the scatter in the TF-relation if late- 
type galaxies are also included in the sample. Because the 
scatter in the TF-relation is often much smaller than the 
indicated values, one may conclude that the TF-relation is 
not simply tracing the connection between luminous mass 
and dynamical stellar mass, and that the amount of dark 
matter is varying systematicly with galaxy color. But in 
short, a shift along one vector in Figs ^ and || is sufficient 
to bring the centers of most T<6 galaxies to one point, 
at least two vectors are needed to bring the centers of the 
T> 6 galaxies to one point. 

A large number of galaxy formation and evolution the- 
ories predict metallicity gradients (for references see e.g. 
VE and Matteucci 1986 , 1992] ) and age gradients (Kenni- 
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cutt 1989 ; Dopita & Ryder 1994 and references therein) 
in galaxies. Several predict both at the same time, like 
viscous galaxy evolution mod els (Lin & Pringle 1987 ; 
Sommer-Larsen & Yoshii 1990 ), the models of Wyse & 



Silk ( 1989] ) and various gas infall and galaxy merger mod- 
els. The young, low metallicity late-type systems observed 
here would be ideal merger candidates to replenish the 
outer regions of large galaxies with nearly unprocessed 
gas from which new stars can be formed. A detailed study 
of all possible models, investigating time scales and metal- 
licity ranges involved, is beyond the scope of this work. 

A key question that is not addressed in this investiga- 
tion is whether the color gradients originate in the arm 
or in the inter-arm region, or are present in both regions. 
High resolution and high signal-to-noise observations are 
needed to solve this question, and the present data set is 
not suited for such an investigation. A detailed study of 
a few large nearby galaxies is in preparation (Beckman et 
al., private communication). 

Another question not addressed in this paper is the 
vertical structure of dust and stellar populations. In a sim- 
ilar approach as followed here, color maps of edge-on gal- 
axies are ve ry well suit ed for such studies (e.g. Wainscoat 
199(j ). Unfortunately, so far most 



et al. 1989, Just et al. 



of these studies have not included scattering in the dust 
models, but the models of Just et al. are most promising, 
as they do include a link between the stellar and dynam- 
ical evolution of the galaxies. 

6. Conclusions 

The stellar and the dust content of a large sample of gal- 
axies was investigated using the color profiles of these gal- 
axies. Data in four optical and two near-IR passbands were 
combined simultaneously to derive structural properties of 
the sample as a whole, rather than for individual galaxies. 
The main conclusions are: 

— Almost all spiral galaxies become bluer with increasing 
radius. 

— The colors of galaxies correlate strongly with surface 
brightness, both within and among galaxies. The mor- 
phological type is an additional parameter in this rela- 
tionship, because at the same surface brightness late- 
type galaxies are bluer than early-type galaxies. 

— Realistic 3D radiative transfer modeling indicates that 
reddening due to dust extinction cannot be the major 
cause of the color gradients in face-on galaxies. The 
predicted color vectors in color-color space are not 
compatible with the data, unless the assumed scat- 
tering properties of the dust are entirely wrong. 

— The color gradients in the galaxies are best explained 
by differences in SFH as function of radius, with the 
outer parts of galaxies being on average much younger 
than the central regions. This implies that the stel- 
lar scalelength of galaxies is still growing. The central 
stellar populations in a galaxy must have a range in 



metallicities to explain the red central colors of the 
galaxies. 

— A consequence of the population changes implied by 
the color differences in and among galaxies is that there 
are large changes in M/L values in and among galaxies. 
These changes in M/L make the missing light problem 
in spiral galaxies as derived from rotation curve fitting 
even more severe. 

— The H and K passbands are recommended for stan- 
dard candle work and for studies depending on M/L 
ratios in galaxies. 
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A. Monte Carlo radiative transfer simulations of 
light and dust in exponential disks 

The modeling of dust extinction in extragalactic systems 
has a long history. The effects of scattering are often ig- 
nored or are assumed to be no more than a scaling factor 
(Disney et al. 1989 ; Huizinga 1994 ), which is correct as 
long as one is not looking at wavelength dependent ef- 
fects. Therefore, scattering is included in most studies in- 
vestigating the wavelen gth d ependent effect s of ex tinction 
(e.g. Kylafis & Bahcall]l987|; Bruzual et al. |1988|, Witt et 



al. |1992[ Byun et al. |1994| ) 

The Monte Carlo approach followed here is a some- 
what "brute force" method, but it has the advantage that 
it can be used for luminous and dust geometries contain- 
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ing little symmetry. Even though smooth light and dust 
configurations are applied here, this method can easily be 
extended to include dust clouds and spiral structure. 

A.l. The mathematical method 

A Monte Carlo radiative transfer code was used to cal- 
culate the light and color distribution of disk galaxies as 
seen by a distant observer. The Monte Carlo principle as 
applied to exti nction in gaseous nebulae is described in 
detail by Witt (|l977j ). In the computer model the paths 
of many (~ 10 6 ) photons were followed as they traveled 
through the absorbing and scattering dusty medium. At 
great distance the photons were collected to produce a 
galaxy image which was used for further study. 

The trajectory of a photon through a dusty medium 
is determined by random processes. One can character- 
ize these processes by a probability function p(x) on an 
interval (a, b) such that 



f P (Od£ = l. 

J a 



(Al) 



Using a random number generator, which produces a 
random number R distributed uniformly in the interval 
0< R < 1, one can simulate an event with frequency p{x)dx 
in the interval (x, x + dx) by requiring 



J a 



(A2) 



This equation was used to simulate the birthplace and 
direction of photons, as well as the scattering properties 
of the dust. In the following sections, each R will denote 
a new random number in the interval < R< 1. 

A. 2. The creation of photons 

This section describes the creation of photons and the ran- 
dom processes involved in this creation. Rather than cre- 
ating photons in space with a certain density distribution 
such as in a real galaxy, the photons were instead created 
uniformly in space and given an initial intensity weight to 
produce the exponential light profile. These initial weights 
were reduced by absorption as the photon moved through 
the dusty medium. 

The models consisted of three-dimensional distribu- 
tions of stellar light and dust that were chosen indepen- 
dently. For the stellar light, an axisymmetric disk-like dis- 
tribution was used, with an exponential intensity behavior 
in both radial (r) and vertical (z) directions: 



I(r,z,< 



-(r/h B + \z\/z B ) 



(A3) 



where h s and z s are the scalelength and the scaleheight of 
the stellar distribution respectively, with the coordinate 
system is as defined in Fig. The luminosity was trun- 
cated at seven scalelengths and heights. The exponential 



behavior of the radial light distribution of disks in spiral 
galaxies is well established, but the vertical light distribu- 
tion is more controversial, because extinction effects make 
measurements difficult. I have used the exponential law 
vertically inste ad of (for instance) the sech or sech 2 laws 
(van der Kruit 1988] ) , because the nearly unobscured light 
distribution as obtained with near-IR observations often 
can be well desc ribed by su ch an exponential law (Wain- 
scoat et al. 1989 Aoki et al. 1991 ). It is not expected that 
the results obtained here will change significantly if other 
plausible vertical light distributions are used. Each "pho- 
ton beam" that was created received an initial weight ac- 
cording to the J(r, z, <j>) of Eq. flA3|). 

Photons were created at a certain position in the model 
galaxy using the distribution functions 

r = 7Rh s z = 7(R- 0.5)2z s <p = 2Rir, (A4) 

and thus their initial position in cartesian coordinates was 



r cos(4>) 
r sin(0) 

z 



(A5) 



The distribution functions of Eq. ( |A4| ) create photons uni- 
formly in cylindrical coordinates, but not in cartesian 
space. The density of created photons is much higher near 
the center than in the outer regions and a correction factor 
is needed. To still get an exponential behavior in flux, an 
additional weight factor of 2itr was given to the intensity 
of each created "photon beam" . These distribution func- 
tions were chosen, because the final model images were 
azimuthally averaged just as the real observations and in 
the absence of dust give these distribution functions an 
equal number of photons at each radius in the face-on 
case. 

At creation, the initial flight direction of the photons 
was specified by the functions 



u = (2R- 1)tt 6=(R- 0.5)tt. 

The corresponding directional cosines are 



A3 



^cos(w) cos(6>) 
sin(w) cos(0) 
sin(0) 



(A6) 



(A7) 



Again, since this would not give a uniform distribution of 
flux density in each direction, an extra weight of the form 
cos(#) was introduced. The initial intensity of a "photon 
beam" at creation was thus 



io(r,z,0,w,0) = 2^rcos(#)e- (r/,l3+|z|/Zs) . 



(A8) 
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Fig. 12. The definition of the coordinate system. A trajectory with a scattering is indicated for one step size. A "photon beam" 
starts at x \d in the direction defined by (w, 9), scatters off a dust particle over an angle a, and ends at x ne w 



A. 3. The dust properties 

Let us first define the notation for the basic equations of radiative transfer. A beam of light with intensity / is followed 
while traveling along a straight line through a dusty medium. Since this medium removes a fraction k ds of the incident 
intensity for each distance ds traveled, the full equation is 



dl(s) 
ds 



-k(s)I(s) + J(s), 



(A9) 



where J(s) are the new sources of intensity at point s in the travel direction, k is the extinction coefficient and consists 
of an absorption and a scattering part, and their relative importance is normally expressed by the albedo (a) 



K — fc a "T K>s 
a = k s /k. 



(A10) 
(All) 



While following the light beam through the dusty medium, no new photons were added to the initial beam of photons 
that was created with intensity 1(0), because they are accounted for by the repeating Monte Carlo principle, and 
J(s)=0. For this single "photon beam" the radiative transfer equation is now 



ds 



= -<s)I(s) 



with formal solution 



The integral in Eq. ( A13 ) defines optical thickness 
t —I k(s) ds. 



(A12) 



(A13) 



(A14) 
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Since the beam was followed while it scattered through the medium (so n o longer on a straight line), there were no 
losses by scattering along the path of the beam. One still can use Eq. ( A13 ), but with K a (s) = (1— o)k(s) instead of k(s) 
as long as the line integral is taken along the path traveled. The extinction coefficient and the albedo are wavelength 
dependent, which creates reddening of stars and could cause the observed color gradients in the galaxies. 

To calculate the absorption along the path traveled, the dust distribution must be defined. The dust extinction 
coefficient in the models was taken to have a similar distribution as the stellar light, but with a scalelength and 
scaleheight that were chosen independently: 



K\(r,z,( 



-(r/h d + \z\/z d ) 



(A15) 



where K\(r, z, <fi) is the local extinction coefficient at wavelength A and ko,a is the extinction coefficient at (r — 0, z — 0). 
In this article the optical depth of a system is denoted by the integration of Ky (k in the y-passband) along the 
symmetry axis from (r — 0, z — ~ oo) to (r = 0, z = oo) 



tq,v 



Ko, y e~ |2|/Zd dz = 2 



KQ.VZd- 



(A16) 



Using Eq. ( A13 ) one can calculate that a point source located behind the center of the galaxy will have suffered an 
extinction in the T^-passband for an observer located at its other pole of 



I(z = oo) = J(0)e" 



(A17) 



assuming that all absorbed and scattered photons are lost for a point source. 

Three parameters are needed to describe dust properties, the relative extinction [t\/tv), the albedo (a\) and the 
scattering phase function & x (a), and all are dependent on wavelength (A). As there are almost no direct measurements 
of the dust p roperties in other g alaxie s, I decided to use the (also poorly determined) Galactic values. Studies by Knapen 
et al. ( |l99l| ) and Jansen et al. ( |1994|) seem to indicate that at least the Galactic extinction curve is applicable to some 
other galaxies, the other two dust properties have never been measured in extragalactic systems. 



For relative extinction, the Galactic extinction properties of Rieke & Lebofsky (1985) were used. The values for 
the albedo were drawn from Bruzual et al. (1988), as were the values for the scattering asymmetry parameter g x . This 
asymmetry parameter enters in the scattering phase function suggested by Henyey & Greenstein (1941) 



$ A (cos(a), 5A ) = [(l-^)/47r](l 



9t 



2g x cos(a))- 3 / 2 



(A18) 



where a is the scattering angle between the incident and the deflected photon. The function is such that g x — (cos(a)), 
and thus — l<g x <l. Forward scattering dominates for <?a>0 and <7a = results in isotropic scattering. To determine 
a ran dom deflection angle for a photon with the probability function $a Eq. (A2) can be used and one finds (Witt 
1977[ ) 



a = arccos({(l + g 2 x ) - [(1 - g 2 x )/(l - g x + 2g x R)} 2 } /2g x ). 

The values used for t\/tv, a Xl and g x for the different photometric passbands are listed in Table [l]. 



(A19) 



A. 4- The numerical method 



Each photon beam that was created using Eqs. (A4) and (|Aq ) was stepped through the dusty medium. The step 
size was determined by the local dust density and was chosen in such way that it would give an optical thickness 
(absorption plus scattering) of 0.03 if there was a constant dust density along the step. This means that the chance of 
multiple scattering is of order (0.03a>,) 2 along a step. This is negligible, as it should be, because only single scattering 
was incorporated along each step. It is trivial to show that the step size has to be of order 



As = om/ K (x) = 0m/(K O e- (r / h ' + ^/ z "'>) 

and the new position is (if there is no scattering) 



(A20) 



AxAt 



(A21) 
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For the actual absorption along a step, one has to use Eq. ( A13 ) with k q = k(1 — a) instead of n. The integral was 



at each step approximated by four-point Gauss-Legendre quadrature: 

r x+AxAs 4 



TAs 



/>X-\-£AXL±S ^ 

/ K a (s)ds w 0.5As^Twj k(w + 0.5(1 + Xj)Ax)(l -a), (A22) 



where the values for Wj and Xj can be found in any handbook on numerical analysis. When there was no scattering 
during the step, the n(x4) calculated in Eq. ( |A22| ) was used to determine the step size with Eq. ( |A2C| ) in the next step. 

If scattering did occur, the step proceeded slightly differently. The chance of scattering during a step can be 
simulated by requiring 

(A23) 



In that case the directional cosines of Eq. (A7) were changed by an angle given by Eq. ( A19). This left another angle 



of freedom for the azimuthal change of the direction which was generated by 

/3 = 2nR. (A24) 

Rotating the initial directional angles uj,9 over the scattering angles a, (3 gives the new directional cosines 

^cos(w) cos(#) cos(/3) +(sin(w) cos(a) + cos(u;) sin(0) sin(a)) sin(/3) 
Ax now = | sin(w) cos(0) cos(/3) +(— cos(w) cos(a) + sin(w) sin(0) sin(a)) sin(/3) | . (A25) 
sin(0) cos(/3) — cos(9) sin (a) sin(/3) 



When Eq. (A23) was satisfied and scattering occurred in a step, the photon proceeded a fraction / = R of step size 



As in the old direction, before following the new direction 

Xncw = Sold + (/Ax id + (1 - f)Ax uew )As. (A26) 
A. 5. Projection on the sky 

The steps in the previous paragraph were repeated until |x n cw| > 10h s . All ta s of the different steps were added, 
and the final intensity of the photon beam was Tend = ioe -Ttot ■ The final intensities were projected on the sky as if 
the galaxy was being observed from infinity. Due to the axisymmetric nature of the models, one can ignore the tu 
dependence of the exit direction and rotate all photons over angle lo as if they leave the galaxy in the same direction. 
The projection on the (yi, 2/2)-plane becomes 

-_( -sinH cosH 0\ 



■ sin(0) cos(cj) — sm(ff) sin(oj) cos(^) 

This procedure of creating photons, stepping through the medium and projecting the exiting photons on the sky was 
repeated for two million photons. The photons were binned into pixel images in (2/1,2/2) direction and the binning 
in viewing angle 9 direction resulted in ten model images from edge-on to face-on. The binning of photons in the 9 
direction was in equal steps of b/a, where b/a is the axial ratio of an inclined circle. This has the advantage that each 
image has approximately equal flux, at least in the case of no extinction. For the analyses discussed in this paper 
only the face-on images were used, which contained the photons that exited with an angle between 90° and 71.8° (i.e. 
l>6/a>0.95). 

A. 6. Testing and the results 

The computer program was tested extensively to ensure that the results were reliable. First, a point source emitting all 



its photons to a pole was placed at the center of the galaxy and the exiting flux was correctly described by Eq. ( AI7 ). 
Tests showed that in the case of no absorption all images at different viewing angle had an equal flux. Tests were made 
to ensure that the scattering angle phase function was produced correctly and that the angle between the incident 
and the deflected photon was indeed the required angle. One of the most important tests was the independence of the 
number of scatterings on step size. If the step size had been too large or if the calculation ta s had been incorrect this 
could not have been the case. The model profiles were tested against the analytic profiles of DDP for the cases with 
and without scattering. In the case of no scattering (all scattered photons were removed to mimic absorption) the 
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agreement was perfect. Figure [| shows an excellent agreement between the two models even with scattering, as long 
as the extinction is high. All test results were correct to within the statistical noise. 

For the extraction of the (color)profiles, the same programs as in Paper I were used and the model results can be 
compared directly with the data. Once a T y and a set of scalelengths and heights have been chosen for a particular 
model, Ko \ is determined by Eq.( A16) and the corresponding optical depths for the other passbands are determined 
by Table [j]. This model is also useful to study inclination-dependent effects of extinction. 
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